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INI KOOUCT ION  AND  SUMMARY 


The  basic  requirement  of  feedback  systems  Is  to  achieve  certain  desired 
levels  of  performance  and  also  to  be  tolerant  of  uncertainties.  Performance 
levels  concern  such  things  as  command  following,  disturbance  rejection,  sen¬ 
sitivity,  etc.,  while  uncertainty  tolerances  deal  with  the  Inevitable  dif¬ 
ferences  which  exist  between  a  physical  plant  and  Its  mathematical  design/ 
analysis  model.  As  discussed  In  various  textbooks  and  references  e.g., 
[1,2],  these  two  aspects  of  the  feedback  problem  lead  to  fundamental  trade¬ 
offs  and  compromises  which  motivate  the  entire  body  of  feedback  theory. 

An  essential  difficulty  In  the  theory  has  been  to  capture  both  the  perfor¬ 
mance  and  uncertainty  aspects  of  feedback  In  a  single  problem  statement. 

Thus  we  have  optimization  theories  which  emphasize  performance,  robustness 
theories  which  emphasize  uncertainties,  and  a  host  of  ad  hoc  tools  which 
attempt  to  compromise  the  two. 

In  this  report,  we  propose  a  problem  formulation  which  captures  both  aspects 
of  feedback  under  the  umbrella  of  what  we  will  call  the  "block-diagonal 
bounded  perturbation  (BDBP)  problem."  The  solution  to  this  problem.  Intro¬ 
duced  In  [3],  Involves  a  generalization  of  the  ordinary  singular  value 
decomposition  (SVD).  It  provides  a  reliable,  nonconservative  measure  to 
determine  whether  both  the  performance  and  robustness  requirements  of  a 
feedback  loop  are  satisfied.  This  measure,  which  we  will  call  the  struc¬ 
tured  singular  value  (SSV)  and  denote  by  the  symbol  u,  serves  as  an  essen¬ 
tial  analysis  tool.  Synthesis  tools  based  on  the  structured  singular  value 
are  under  development. 

The  report  Is  organized  Into  four  major  sections.  Section  2  formulates  the 
robustness  and  performance  aspects  of  feedback  as  a  block-diagonal  bounded 
perturbation  problem.  This  problem  Is  solved  In  Section  3  using  the  new 
structured  singular  value  concept.  Performance  Implications  of  these 
results  are  then  examined  In  Section  4,  which  Includes  a  fundamental  theorem 
relating  performance  In  the  face  of  uncertainty  to  the  SSV,  as  well  as  an 
example.  Section  5  contains  concluding  comments. 


FEEDBACK  ANALYSIS  AS  A  BLOCK-OIAGONAL  BOUNDED 
PERIURBATION  PROBLEM 


This  section  formulates  the  basic  feedback  problem  of  achieving  performance 
in  the  face  of  uncertainties  as  a  stability  problem  in  the  presence  of  block- 
diagonal  bounded  perturbations.  The  formulation  Involves  cone-bounded 
transfer  functions  as  basic  building  blocks,  in  terms  of  which  both  the 
robustness  and  performance  aspects  of  feedback  can  be  characterized. 

Basic  Building  Blocks:  Cone  Bounded  Transfer  Functions 
Throughout  this  report,  we  will  deal  with  multi-variable  feedback  systems 
whose  models  are  linear,  time  invariant,  and  finite  dimensional.  Hence, 
they  can  be  represented  by  transfer  function  matrices  with  rational  ele¬ 
ments.  The  robustness  and  performance  properties  of  these  systems  will  be 
expressed  in  terms  of  a  collection  of  transfer  matrices,  A^s),  1=1,2, 

•••,  m,  which  each  satisfy 

A1  (s)  -  L"1  (s)  eCs)!^  (s)  (1) 

where  L^s)  and  R ^ ( s )  are  constant  transfer  matrices  and  0(s)  is  any 
stable  transfer  matrix  from  a  set  satisfying 

a[9  (Ju)]  <1 

We  will  also  require  that  L^  and  R^  have  no  poles  or  zeros  in  the  open 
right  half  plane.  These  assumptions  assure  that  A^  has  no  rhp  poles. 

Reasons  for  this  restriction  are  discussed  later.  Note  that  the  functions 
which  satisfy  (1)  belong  to  conic  sectors,  as  initially  defined  by  7ames 
([A],  [5])  and  generalized  by  Safonov  ([6],  [1]).  Their  sector  centers  are 
zero,  and  their  sector  radii  are  characterized  by  L(s)  and  R(s).  We  will 
use  such  cone  bounded  transfer  functions  as  basic  building  blocks  In  a  com¬ 
bined  robustness/performance  characterization  of  feedback  systems. 

Robustness  Characterization 

The  use  of  cone  bounded  transfer  functions  to  characterize  robustness  has 
been  a  central  theme  in  many  recent  references,  including  [8]  where  such 
transfer  functions  were  Inserted  at  the  Inputs  or  outputs  of  a  plant  model 
in  order  to  represent  so  called  unstructured  uncertainties  (modelling  errors 
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with  no  assumed  structure  except  for  known  magnitude  bounds  on  their  trans¬ 
fer  functions).  Necessary  and  sufficient  conditions  were  then  derived  for 
stability  robustness  In  the  face  of  such  uncertainties.  For  example,  a 
stable  feedback  loop  with  plant  G(s)  and  compensator  K(s)  will  remain  stable 
In  the  face  of  all  possible  perturbed  plants  G'(s)  =  [I  ♦  A^]G(s),  with  A^ 
given  by  (1),  If  and  only  If 

5[R^  GK( I  +  GK)-^  ]  <  H  a  >  0  (2) 

Note  that  with  R  and  l  specified,  this  Inequality  Imposes  conditions 

'  -1 
on  the  shape  of  the  closed  loop  frequency  response,  GK(I  +  GK)  ,  which 

must  be  satisfied  In  order  to  assure  robust  stability.  These  conditions 
are  unique  to  the  assumed  form  of  plant  perturbations  (e.g.,  G'  =  (I  +  A)G 
In  the  present  case).  Each  such  assumed  form  corresponds  to  a  specific 
location  where  A  Is  Inserted  In  the  nominal  feedback  loop.  The  location 
for  our  present  case  Is  shown  In  Row  1  of  Table  1.  Other  locations  corres¬ 
pond  to  other  assumed  forms  for  G'  and  produce  different  necessary  and  suf¬ 
ficient  stability  robustness  conditions.  A  representative  set  of  possibili¬ 
ties  Is  summarized  In  the  remaining  rows  of  Table  1.  (Host  of  these  cases 
can  be  found  In  [9]) . 

Table  1  also  Indicates  representative  types  of  physical  uncertainties  which 
can  be  usefully  represented  by  cone  bounded  perturbations  Inserted  at  the 
Indicated  locations.  For  example,  the  representation  G'  =  (I  *  A)G  In 
Row  1  Is  useful  for  output  errors  at  high  frequencies,  covering  such  things 
as  unmodelled  high  frequency  dynamics  of  sensors  or  plant.  Including  diffu¬ 
sion  processes,  transport  lags,  electro-mechanical  resonances,  etc.  The 
representation  G*  »  G(I  ♦  A)  In  Row  2  covers  similar  types  of  errors 
occusrlng  at  the  Inputs.  Both  cases  should  be  contrasted  with  Rows  4  and  5 
which  treat  G'  =  (I  ♦  A)_1G  and  G'  *  G(I  ♦  A)-1.  These  representations  are 
more  useful  for  variations  In  modelled  dynamics,  such  as  low  frequency 
errors  produced  by  parameter  variations  with  operating  conditions,  with 
aging,  or  across  production  copies  of  the  same  plant.  Discussion  of  still 
other  cases  Is  left  to  the  table.  Note  from  the  table  that  the  stability 
requirements  on  A  do  not  limit  our  ability  to  represent  variations  In 
either  the  number  or  location  of  rhp  singularities. 
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The  most  significant  thing  to  understand  about  Table  1  Is  that  the  stability 
robustness  conditions  shown  are  sufficient  to  assure  stability  only  If  all 
the  uncertainties  occur  at  the  Indicated  locations  and  none  occur  else¬ 
where.  In  order  to  use  the  conditions  directly,  therefore,  designers  are 
obliged  to  reflect  all  known  sources  of  uncertainty  from  their  known  point 
of  occurance  to  a  single  reference  location  In  the  loop.  Such  reflected 
uncertainties  Invariably  have  a  great  deal  of  structure  which  must  then  be 
"covered  up"  with  a  larger,  arbitrarily  more  conservative  perturbation  In 
order  to  maintain  a  simple  cone  bounded  representation  at  the  reference 
location.* 

Alternatively,  designers  could  choose  to  treat  uncertainties  occurring  at 
several  different  locations  In  the  feedback  loop  as  a  single  uncertainty 
occurring  at  one  location  In  a  larger  feedback  loop.  To  be  specific  about 
this  alternative,  let  A^ ,  1=1,  2,  •••,  m,  denote  a  collection  of  such 
uncertainties  positioned  at  location  8,^,  1=1,  2,  •••,  m.  Note  that  at  each 
the  feedback  loop  has  an  Input,  where  It  receives  the  signals  from  A^ , 
and  also  an  output,  where  It  supplies  signals  to  A^ .  Let  be  the 
transfer  function  matrix  between  these  two  sets  of  signals.  Further,  let 
denote  the  transfer  matrix  between  the  Inputs  at  location  and 
the  outputs  at  location  1  .  Then  the  block-structured  matrix 

M  =  {  }  (3) 

represents  all  Interactions  of  the  feedback  loop  with  Its  uncertainties,  and 
Indeed,  the  block-diagonal  bounded  perturbation  diagram  In  Figure  1  Is  an 
equivalent  representation  of  the  loop.  Here  we  have  A  =  d1ag(A.j,  A2,***,Am). 

Note  that  the  feedback  elements  In  this  larger  loop  are  zero  In  the  absence 
of  uncertainties.  Hence,  M  will  be  a  stable  "plant"  whenever  the  original 
nominal  loop  Is  stable.  As  an  example  of  this  representation,  consider  the 
system  In  Figure  2.  This  system,  with  two  uncertainties  present  simulta¬ 
neously,  the  first  from  Row  2  and  the  second  from  Row  4  of  Table  1,  Is 
described  by  the  following  M  matrix. 

*  By  "arbitrarily  more  conservative,"  we  mean  that  examples  can  be  construc¬ 
ted  where  the  degree  of  conservatism  Is  arbitrarily  large.  Of  course,  other 
examples  exist  where  It  Is  quite  reasonable. 


M 


(4) 


(I  ♦  KG)_1KG  (I  *  KG)"1 

(I  *  GK)-1G  (I  ♦  GK)'1 

Given  the  equivalent  system  In  Figure  1  with  A^'s  characterized  by  equa¬ 
tion  (1),  It  follows  from  the  Small  Gain  Theorem  [10]  that  the  loop  remains 
stable  In  the  presence  of  these  uncertainties  If 

a  [R(j«)M(jG))L(juJ‘1  ]<1  H>0  (5) 

where  R  =*  d1ag(R.  R„  •••  R  )  and  L  =  d1ag(L,  l.  •••  L  ) 
i  i  m  1  c  m 

This  condition  provides  an  alternate  test  for  stability  robustness.  Like 
the  procedure  of  reflecting  all  uncertainties  to  one  reference  location, 
however,  the  new  test  can  be  arbitrarily  more  conservative  because  It 
Ignores  the  known  block-diagonal  structure  of  the  uncertainties  In  Figure  1. 

The  objective  of  our  results  In  this  report  Is  precisely  to  reduce  the  con¬ 
servatism  of  robustness  and  performance  tests  for  block  diagonal  structures 
such  as  Figures  1.  We  do  this  by  Introducing  a  generalized  notion  of  the 
maximum  singular  value  for  block-diagonal  structures.  This  generalization 
Is  developed  In  Section  3.  It  Is  called  the  structured  singular  value  (SSV) 
and  Is  denoted  by  the  symbol  p.  It  yields  the  following  necessary  and 
sufficient  conditions  for  robust  stability  of  the  BDBP  problem: 

u[R(Ju)M(ju)L1  (Jo)]  <  1  V-  o  >  0  (6) 

This  represents  our  extension  of  the  Small  Gain  Theorem  which  we  call  the 
Small  p  Theorem. 

Since  all  simultaneous  uncertainties  can  be  put  Into  block-diagonal  form  by 
merely  constructing  the  associated  matrix  M,  the  SSV  allows  us  to  nonconser- 
vatlvely  analyze  simultaneous  occurences  of  uncertainties  anywhere  In  a 
feedback  system.  The  uncertainties  may  be  cone  bounded  errors  of  Individual 
components  of  the  system  (SISO  or  MIMO),  they  may  be  Individual  parameter 
variations  In  the  model,  or  even  polynomial  approximations  of  parameters 
entering  nonllnearly.  In  fact,  the  only  restrictions  which  remain  Is  that 
all  variations  must  be  allowed  to  be  complex.  Pure  real  variations  or  pure 
Imaginary  variations  cannot  be  separated  Into  Individual  blocks. 
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Performance  Characterization 

The  ability  to  treat  simultaneous,  structured  uncertainties  also  offers, 
almost  as  a  free  byproduct,  the  ability  to  deal  simultaneously  with  the  per¬ 
formance  and  robustness  aspects  of  feedback.  This  Is  made  evident  In  Column 
4  of  Table  1,  where  each  of  the  conditions  Imposed  on  feedback  loop  shapes 
by  perturbation  A^  at  location  2^  Is  given  a  performance  Interpreta¬ 
tion.  For  example,  the  perturbations  In  Row  4  Impose  requirements  (through 
l  and  R)  on  the  shape  of  the  function  (I  +  KG)-1.  This  function  Is,  of 
course,  the  classical  (output)  sensitivity  function  of  the  feedback  loop. 
Small  values  over  some  frequency  range  guarantee  low  closed  loop  sensitivity 
to  open  loop  variations  and  low  command  following  errors  to  output  commands 
over  that  range.  A  particular  specification  on  these  performance  parameters 
can  thus  be  imposed  on  a  design  by  Introducing  a  "fictitious  uncertainty"  at 
the  location  In  Row  4  with  cone  bounds  R  and  L  selected  to  meet  the  perfor¬ 
mance  requirement. 

To  Illustrate  how  such  fictitious  uncertainties  actually  enforce  performance 
specs,  consider  the  simple  case  where  a  single  true  uncertainty,  say  Ar 
from  Row  2,  and  a  single  fictitious  (performance)  uncertainty,  cay  Ap 
from  Row  4,  are  specified  for  our  feedback  system.  Let  the  structured  sin¬ 
gular  value  condition  (6)  be  satisfied  for  the  corresponding  M  matrix 
(equation  4).  Then  the  system  remains  stable  In  the  face  of  Ar  and  Ap 
occurring  simultaneously.  Obviously,  It  will  also  remain  stable  for  Ap 
with  Ar  *  0.  This  means  that  the  nominal  system  must  satisfy  the  per¬ 
formance  condition 

o[Rp(I  +  KG)'1  Ip1  ]  <  1  *  o  >  0  ,  (7) 

because  the  latter  Is  also  a  necessary  and  sufficient  condition  for  robust 
stability  with  Ap  only.  This  much  Is  straightforward.  What  Is  not  so 
evident  but  much  more  Important  Is  that  Condition  (7)  Is  also  satisfied  for 
all  perturbed  feedback  loops.  That  Is,  for  all  true  plants  G'=G(I*Ar) 
we  have 

o  [  Rp  ( I  ♦  KG')1  Ip"1]  <  1  -V  a  >  0.  (8) 

Hence,  the  performance  spec  Is  satisfied  In  the  face  of  all  possible  true 
uncertainties.  A  proof  of  this  consequence  of  the  structured  singular  value 
condition  Is  left  to  Section  4. 


SECTION  3 

THE  STRUCTURED  SINGULAR  VALUE  AND  THE  SMALL  p  THEOREM 

We  have  discussed  how  the  problem  of  analyzing  performance  In  the  face  of 
structured  uncertainty  can  be  expressed  as  a  BDBP  problem.  We  noted  that 
standard  singular  value  tests  applied  to  the  BDBP  can  be  excessively  conser¬ 
vative  because  they  Ignore  the  block  diagonal  structure.  A  more  general 
non-conservative  test  (the  Small  p  Theorem)  Is  developed  In  this  section 
which  removes  this  limitation.  By  non-conservative  we  mean  providing  a 
necessary  and  sufficient  condition.  The  test  Is  expressed  In  terms  of  a  new 
measure,  the  structured  singular  value  p.  This  section  begins  with  review 
of  the  results  In  [3],  where  p  was  Introduced. 

To  provide  a  more  precise  description  of  block  diagonal  perturbations,  let 
K  =  (m1  ,m2,***,mn,  k1  ,k2,*»*,kn)  be  a  2n-tuple  of  positive 
Integers.  All  the  definitions  that  follow  depend  on  K,  but  to  simplify  not¬ 
ation  this  dependency  will  not  be  explicitly  represented.  Let 


k  =  £  m,  k,  and  m  =  £ 

J=1  J  J  J=1 


£  m,  . 
J=i  J 


For  each-6c[0,oo),  let  X  efH(k)  be  a  set  of  6-norm-bounded 

6 

block-diagonal  matrices  defined  by 


xw=  {d1ag(A  A  Y  &2,& 2,***,A2,  a3,  An  l,  A  rf ***,fln) 

|Aje7fl(kj)  and  5(Aj  )  <6  for  each  j  =  l,2,***n} 


Let  X  =  U  X,  be  the  set  of  all  such  matrices  with  no  restriction  on  the  norm, 

*  4.1  1 

'LL  be  the  set  of  block  diagonal  unitary  matrices, 

U=  U(k)/|X1 

and  s>  the  set  of  real  diagonal  matrices  such  that 
»  *  I  “'*><■',  V  Vk-  —  \  \  dm,.  1  'k,-  li  \  >  (10) 

d^c  R+*  (0,  <*») } 
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IS  BLANK 


(11) 


What  Is  desired  Is  a  function  (depending  on  “X) 

y.:Tn(k)  -»  [0,  ®) 

with  the  property  that  for  V  Me7fl.(k) 

det(  I*MA)  4  0  ¥■  AeX& 
Iff  8u(M)£l 


(12) 


This  could  be  taken  as  a  definition  of  y..  Alternatively,  y.  could  be 
defined  as 

(0  If  no  AeX  solves  det(I+MA)  =  0 

(13) 

IaIx  (0(A)  I det( I+MA)  =  0}|  _1  otherwise 

This  definition  shows  that  a  well-defined  function  satisfies  (12).  It  pro¬ 
bably  has  little  additional  value  since  the  optimization  problem  Involved 
does  not  appear  to  have  useful  properties. 


Using  these  definitions,  the  following  useful  properties  of  y.  are  easily 
proven. 


1)  u(aM)  =  |  a  |  y.(M) 

2)  u(I)  =  1 

3)  u(AB)  £  5(A)u(8) 

4)  y.( A)  *  o( A) 

5)  If  n=l  and  m.^1  then  v(M)  «  ®(M) 

6)  If  n=.1 ,  k  «1,  then  k=m^ , 

X8  =  {\Il\cf|\|  <  6}  and  y(M)  =  p(M) 

7)  If  AeX  ,  UeOt  then  UAeXt  and  AUeX. 

0  o  d 

8)  For  ¥AeX  and  ¥-0eK 

CO 

9)  For  -V-  UeU  and  Me-m(k) 


¥  Metn(k) 

¥  AeX,  Be7»l(k) 
¥  AeX& 

¥  Me7)?(k) 

¥•  Mem(k) 
DAD_1=A 


U(HU)  =  u(UM)  *  y(H) 

10)  For  ¥-  Dc B  and  Mem(k) 

y.( DMO^)  *  u(H) 

11)  max  p(UM)  £  u(H)  £  Inf  a  (OMD'1  )  for -V  Me?n.(k) 

UelC  De# 
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Properties  5)  and  6)  show  that  the  structured  singular  value  has  as  special 
cases  both  the  spectral  radius  and  the  maximum  singular  value.  Property  9) 
means  that  u  Is  ^.-Invariant. 

The  most  Important  results  from  [3]  are  the  following,  which  deal  with  the 
bounds  In  property  11): 

a)  The  left-hand-side  Inequality  In  (11)  Is  always  an  equality. 

This  expresses  y.  In  familiar  linear  algebraic  terms,  but  the 
optimization  problem  Involved  may  have  multiple  local  maxima. 

b)  The  right-hand-side  Inequality  In  (11)  Is  an  equality  when 
there  are  three  or  fewer  blocks,  and  the  blocks  are  not 
repeated.  The  blocks  themselves,  and  therefore  M,  may  be  of 
arbitrarily  large  dimension.  A  tedious  but  straightforward 
computation  shows  that  the  optimization  problem  Involved  Is 
always  convex  [11].  Furthermore,  the  minimization  Is  over 
only  n-1  parameters  for  n  blocks,  Independent  of  block  size, 
making  this  an  attractive  alternative  to  a). 

Note  that  the  transformation  DM0-1  Is  simply  a  rescaling  of  the  Inputs  and 
outputs  of  M.  The  SSV  Is  Invariant  with  respect  to  such  rescaling  (property 
10),  while  singular  values  do,  of  course,  vary  with  rescaling.  This 
Implies,  for  example,  that  the  ad  hoc  method  of  performing  a  change  of  units 
can  reduce  the  conservatism  associated  with  singular  values.  For  some  time 
we  have  been  using  Osborne's  technique  [12],  which  minimizes  the  Frobenlus 
norm  of  DMD  1  to  compute  frequency-dependent  D  matrices.  We  now  have  new 
algorithms  which  compute  0  to  directly  minimize  5(DMD~  ). 

Numerical  software  for  computing  y.  has  been  developed  using  algorithms 
based  on  these  results.  In  addition  to  using  this  software  to  analyze  some 
simple  feedback  designs,  test  runs  have  been  made  on  a  large  number  of 
psuedo-random  matrices.  It  appears  that  the  global  maximum  In  a)  Is  often 
easily  found,  although  a  simple  gradient  search  Is  Inadequate.  Also,  the 
bound  obtained  In  b)  appears  to  be  quite  good  (to  within  15%)  for  cases  of 
more  than  3  blocks.  These  observations  are  most  encouraging,  especially 
considering  the  experimental  and  preliminary  nature  of  the  software. 


There  are  essentially  two  direct  applications  of  singular  values  to  the  BDBP 
problem,  which  provide  bounds  for  u: 

1)  Ignore  the  block  diagonal  structure  and  compute  o(M).  This 
gives  an  upper  bound  for  u- 

2)  Treat  each  perturbation  one  at  a  time.  Compute  the  largest 
maximum  singular  value  for  each  of  the  corresponding  diagonal 
blocks.  This  gives  a  lower  bound  for  u- 

The  gap  between  these  two  bounds  may  be  arbitrarily  large. 

An  extension  to  1)  was  proposed  by  Lehtomakl  ([9], [13]),  who  uses  the  sin¬ 
gular  vectors  for  <5(M)  to  sharpen  the  bound.  Lehtomakl' s  method  checks 
for  structure  but  not  In  the  BDBP  form.  The  optimism  of  2)  can  be  reduced 
by  using  a  method  suggested  by  Freudenberg,  et  al  [14],  who  evaluate  the 
differential  sensitivity  of  the  singular  values  at  one  point  with  respect  to 
perturbations  at  another.  Although  this  method  does  not  apply  to  simulta¬ 
neous,  large  perturbations.  It  can  be  quite  useful  In  Indicating  when  the 
lower  bound  for  u  obtained  by  method  2)  Is  optimistic.  It  should  be  men¬ 
tioned  that  Lehtomakl  and  Freudenberg  did  not  present  their  techniques  In 
the  context  of  the  BDBP  problem. 

The  preceding  discussion  of  p.  and  the  BDBP  problem  has  dealt  with  deter¬ 
mining  the  size  of  the  minimum  structured  perturbation  A  that  causes  I  * 

HA  to  be  nonsingular.  We  are  Interested  In  using  the  structured  singular 
value  to  answer  robustness,  sensitivity,  and  performance  questions  for  mul¬ 
tivariable  feedback  systems.  The  connection  between  y.  and  these  essential 
feedback  properties  Is  provided  by  the  Small  p  Theorem,  which  characte¬ 
rizes  the  stability  robustness  properties  of  a  feedback  system  with  respect 
to  block  diagonal  perturbations.  In  order  to  state  the  Small  p  Theorem  we 
need  the  following  additional  definitions  (all  depending  on  K): 

Let  L,  R  eft.  /]{X  x  t  }  be  such  that  L  and  R  have  no  poles  or 
00 

zeros  In  the  open  right-half-plane.  Then  let 

w  »  tecR/KX^x  e)  |  o(0(S))  £  1  -v  Re( s)  >  0] 


(16) 


and  X  *  {L_1©R|  6e>}. 


(17) 


For  the  BDBP  problem  In  Figure  1,  Me^,(k),  x  Is  the  set  of  allowable 
block  diagonal  perturbations,  and  L  and  R  are  the  weightings  for  the  A 
such  that  3(LAR~  )  £.  1.  We  will  say  the  canonical  system  In  Figure  1 
Is  stable  Iff  I+HA  Is  nonsingular  In  the  closed  right-half-plane. 

Although  this  definition  does  not  distinguish  between  111-posedness  and 
Instability,  It  Is  adequate  for  our  purposes.  We  can  now  state  and  prove 
the  following: 

Theorem  (Small  u) :  The  canonical  system  Is  closed  loop  stable  for  all 
Ac  ^  Iff 

Uc  *  sup  vl(RML_1  )  <  1  (18) 

u 

Proof:  To  prove  the  If  part,  suppose  v.c  <.  1  and  let  A  eX.  Then 
using  Properties  3)  and  11)  and  the  definition  of  X 

sup  p( MA)  *  sup  p(MA)  4  sup  ti(RML-1)  =  u  *  1 . 

Res>0  s*)«  s*Jq  1 

Thus  I  +  MA  Is  nonsingular  for  all  Res  >  0.  Since  A  was  arbitrary, 
the  canonical  system  Is  stable  for  all  A  c 

Conversely,  suppose  y.(RML_1)|  £  1  (u  may  be  ®).  Then  3  0cX,,3, 

&>=Ci>  0  s 

0 

det( I  ♦  RML-18)  %  =  0.  Thus,  3  A  eX  *  det(I  +  MA)  |  =0  and  the 

0)3  Cl)  *  * 

0  0 

canonical  system  Is  not  stable  for  all  A  c  H.  □ 

This  theorem  guarantees  that  If  ii(RMl_1)  Is  less  than  1  at  every  fre¬ 
quency,  then  the  closed-loop  system  Is  stable  for  all  structured  perturba¬ 
tions  Acj£.  Conversely,  If  v(RML-1)  Is  greater  than  or  equal  to  1  at 
some  frequency,  then  there  exists  a  structured  perturbation  AeXthat 
results  In  closed  loop  Instability.  Note  that  a  destabilizing  A  can  be 
expressed  as  L8R-1  for  some  constant  0. 
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As  noted  In  Section  2.  the  Small  y  Theorem  can  also  guarantee  a  pre¬ 
specified  performance  level  by  Including  a  performance  block  In  the  BOBP 
problem.  Furthermore,  this  performance  level  Is  guaranteed  for  all  struc 
tured  perturbations  Aejt.  These  claims  are  made  precise  In  the  next 
section  by  a  corollary  to  the  Small  y  Theorem  that  treats  performance. 


SECTION  4 

PERFORMANCE  IMPLICATIONS 


Suppose  that  the  plant  perturbations  are  given  by 

\  ■  <11a9lAr  V  V  c3Er 

with  corresponding  weighting  matrices  Lr,  Rr,  Interconnection  matrix 

M  ,  and  2n-tup1e  K  =  (m,,  •  •«,  m  ,  k,,  •••,  k  ).  Suppose  that  a  per- 
r  r  1  n  1  n 

formance  specification  Is  given  as 

^RpM^JLp1  )  1  -V  g>  V-Are  Xr  (19) 

Here  M'  Is  a  k  xk  performance  matrix  which  we  desire  to  be  small  (as 

P  P  P  ,  , 

weighted  by  Rp  and  Lp).  Examples  Include  M'p=(I+G'K)  ,  Mp=G’K(I  *  G'K)  , 

etc.,  as  discussed  In  Section  2.  Note  that  M'  depends  on  the  perturba¬ 
tion  Ar,  Indicating  that  this  performance  should  be  met  for  all 
uncertainties . 


Let  M  and  M  denote  the  transfer  function  matrices  between  performance 
pr  rp 

outputs  and  perturbation  outputs  and  between  perturbation  Inputs  and  performance 
Inputs,  respectively.  In  terms  of  these  matrices,  It  can  be  shown  that 


"W  "  Hp  ♦  VV1  +  Hr\> 


-1 


rp 


where  Mp  Is  the  performance  matrix  In  the  absence  of  uncertainties. 
Define 

Kt  -  (1,  nuj ,  •••,  mn,  kp,  ky  •••,  kn) 

Lt  -  d1ag(Lp,  Lr) 

Rt  »  d1ag(Rp,  Rf) 

)+r  'VKJ 
vp  -urn.  kr)) 

>ij  *“H(K»j)  “"Rp  x‘Wr 
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We  have  noted  the  dependence  on  K  here  to  avoid  confusion.  Of  course,  the 
nominal  Interconnection  matrix  M  Is  assumed  to  be  stable. 

For  these  definitions,  the  following  relationship  exists  between  perfor¬ 
mance,  stability  robustness,  and  the  SSV. 

Theorem:  (Robust  Performance) 

Mp'(Ar)  stable  and  a(R  ^  ( A  ^  1  V-  u  and  «  A  r  c  2r  (20 

Iff  )il  ¥  u 

Proof:  It  follows  from  the  Small  y.  Theorem  that 
y(RTHrLT'1)  <  1  ¥  ca 

Iff  |I  ♦  RrNrlT’<le|  >  0,  -V  Res  >  0,  -V  0  c>T 

Iff  ¥  Res  >0,  |I  ♦  R  H  L  9  |  >  0,  ¥  0„  el* 

r  r  r  r  r  r 

>nd  |I  *Rpy  ( ^  e  ^  ^  L -1  0p  I  >  0.  V0peWp 

Iff  ¥  Res  >  0,  |I  *  MfA  r  I  >0  and 

11  *RpV  «4^Li'9p  l>°  V-0pcWp 

Iff  Mp|  (/y)  stable  and  5(RJ*pl  (\)C]p  )  *  1 
for  ¥  u  and  ¥  Ay.  e^.  □ 

We  note  that  this  theorem  extends  the  Small  y.  Theorem's  robust  stability 
results  to  a  composite,  simultaneous  result  on  robust  stability  and  perfor¬ 
mance.  Thus,  given  an  uncertain  plant  model  with  structured  perturbations 
and  a  performance  specification,  we  have  a  necessary  and  sufficient  condi¬ 
tion  In  terms  of  y.  for  satisfaction  of  the  performance  spec  In  the  face  of 
the  uncertainty.  If  the  condition  y  <.  1  Is  met,  then  the  desired  perfor¬ 
mance  Is  achieved  for  all  perturbed  plants.  If  y>  1,  then  there  exists  a 
structured  perturbation  which  causes  the  performance  spec  to  be  violated. 
The  robust  performance  condition  may  be  thought  of  as  arising  from  an  equi¬ 
valent  "fictitious  uncertainty,"  although  this  Interpretation  Is  not 
necessary. 


Example 

In  this  section  the  Ideas  we  have  presented  will  be  Illustrated  by  a  simple 
"textbook*  example.  An  example  similar  to  this  was  discussed  In  [8]. 


The  configuration  for  this  example  Is  shown  In  Figure  2  and  was  discussed 
briefly  In  Section  2.  The  plant  transfer  function  Is 


the  uncertainty  weightings  are  L^L^I,  R^^I,  and  R?  =  r2I,  where 
r^ ( s )  =1  *  s  and  r2(s)  =  1  «•  ^ 


(22) 


Interpreted  In  terms  of  uncertainty  levels,  these  weightings  mean  that 
A^  Is  large  at  high  frequency  and  A2  Is  large  at  low  frequency.  A  design 
that  yields  y.  *  a.  ■¥■  o  has  the  property  that  the  closed  loop  system 
would  remain  stable  for  all  simultaneous  perturbations  such  that 


Larger  perturbations  would  destabilize  the  closed  loop  system.  An  alterna¬ 
tive  Interpretation  In  terms  of  performance  would  be  that  low  the  closed 
loop  system  has  the  following  output  sensitivity 


a[(l  ♦  KG')-1] 

* 


In  the  face  of  all  Input  perturbations  which  satisfy 


otA  ^  —  -V  « 


The  weighting  r^  emphasizes  sensitivity  at  low  frequency. 


(24) 


(25) 


This  example  Is  not  motivated  by  a  physical  design  problem,  and  either 
Interpretation  Is  possible.  We  will  simply  compute  u  for  3  different 
designs  and  compare  them  with  each  other  rather  than  with  a  performance/ 
robustness  specification.  We  will  also  compute  some  singular  value  bounds 
and  u  for  the  case  of  r^  »  r2  «  1 .  No  claims  are  made  about  the 
quality  of  these  designs.  They  are  Intended  merely  to  be  Illustrative. 
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The  M  matrix  for  this  example  Is  given  In  equation  (4)  and 

-1 


U(RML_1)  =  a 


rlTl 

dr2S2G 


d  r^G 


r2S2 


where  T1  =>  (I  *  KG)_1KG  and  »  (I  ♦  GK)-1.  The  compensator 
matrices  for  the  designs  are 

-8  9 


1) 


2) 


3) 


9  -10 


K2  -  0^2 


C3 


where  Q 


10(  s»1 ) 


V 


3s{ s*l 6) 
0 

3/4 
2/3 


9il6s+l) 

32s(s+l) 


-2/3 

3/4 


.  Q, 


3/4  2/3 

2/3  -3/4 


(26) 


(27) 


In  the  first  design,  basically  tnverts  the  plant,  making  KG=GK*7/s  I. 

The  unweighted  singular  values  for  the  two  uncertainties  considered 
Individually  (a(T.j)  and  a(S2))  are  plotted  as  the  lower  two  curves  In 
Figure  3.  The  maximum  of  these  give  a  lower  bound  for  the  unweighted  ii(M). 
The  singular  values  of  M  are  the  upper  two  curves  In  Figure  3.  The  maximum 
of  these  give  an  upper  bound  for  v-  The  unweighted  SSV  Itself  Is  shown 
as  the  plots  labelled  with  a  1  In  Figure  4. 


This  example  Illustrates  that  v  can  equal  the  upper  bound  o(H)  or  lower 
bounds  a(T.j)  and  o(S2)  or  be  anywhere  In  between.  Note  that  this  design 
can  tolerate  either  or  separately  as  large  as  1,  but  simultaneous 
variations  of  less  than  0.1  can  be  destabilizing.  Interpreted  In  terms  of 
performance  for  the  weightings  In  (22),  we  can  see  that  the  design  has  good 
output  sensitivity  (o(r2S2)  ±  1)  for  the  unperturbed  plant.  Small  Input 
perturbations,  however,  can  lead  to  very  poor  output  sensitivity  since 
h(RNL_1)  >>  1  Implies  that  a(r2S'2)  >>  1  for  some  perturbed  S2.  Designs  2  and 
3  have  substantially  Improved  unweighted  SSV's,  as  seen  by  the  plots 
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labelled  2  and  3  In  Figure  4.  It  appears  from  these  plots  that  Designs  2 
and  3  are  uniformly  better  than  Design  1,  with  Design  2  best  In  the 
midfrequency  range  and  Design  3  best  at  the  high  and  low  frequency 
extremes.  The  significance  of  these  differences,  however,  can  only  be 
Interpreted  against  the  given  performance  and/or  robustness  specs.  This  Is 
done  with  the  corresponding  weighted  SSV's  In  Figure  5.  We  see  that 

18  for  Design  1 

a  =  max  u[RHL_1  ]  =  18-5  for  Design  2 

6) 

5.5  for  Design  3 

Thus,  Design  3  offers  the  smallest  a-value  for  equations  (23)  -  (25)  and 
Is  the  best  system  when  judged  against  the  spec. 

This  example  has  Illustrated  several  Important  points: 

1)  Unsealed  singular  value  bounds  provide,  In  general,  poor  estimates 
for  v. 

2)  Stablllty/performance  evaluations  with  perturbations  one  at  a  time 
can  be  highly  optimistic. 

3)  Designs  can  be  meaningfully  compared  only  with  respect  to  some  per 
formance  and  robustness  specification. 


SECTION  5 
CONCLUSION 


This  report  Introduced  an  analysis  technique  based  on  the  Structured  Sin¬ 
gular  Value  u  for  linear  feedback  systems  that  provide  a  reliable,  noncon¬ 
servative  measure  of  performance  In  the  face  of  structured  uncertainty.  The 
Small  y.  Theorem  gives  a  necessary  and  sufficient  condition  In  terms  of  y 
for  stability  of  a  linear  system  with  multiple,  simultaneous,  norm-bounded 
perturbations  of  fixed  but  arbitrary  structure.  The  Robust  Performance 
Theorem  provides  a  similar  condition  for  the  satisfaction  of  performance 
specifications  In  the  presence  of  structured  perturbations.  Some  simple 
feedback  designs  were  presented  to  Illustrate  the  theory. 
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Figure  1.  Feedback  loop  as  a  BD8P  Problem 
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Figure  5.  Weighted  v  for  three  designs 


